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Abstract 

A wide variety of dissipative and fluctuation problems involving a quantum system in a heat 
bath can be described by the independent-oscillator (10) model Hamiltonian. Using Heisenberg 
equations of motion, this leads to a generalized quantum Langevin equation (QLE) for the 
quantum system involving two quantities which encapsulate the properties of the heat bath. 
Applications include: atomic energy shifts in a blackbody radiation heat bath; solution of the 
problem of runaway solutions in QED; electrical circuits (resistively shunted Josephson barrier, 
microscopic tunnel junction, etc.); conductivity calculations (since the QLE gives a natural 
separation between dissipative and fluctuation forces); dissipative quantum tunneling; noise 
effects in gravitational wave detectors; anomalous diffusion; strongly driven quantum systems; 
decoherence phenomena; analysis of Unruh radiation and entropy for a dissipative system. 
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I. INTRODUCTION 



Noise is due to rapid fluctuations in the average of physical quantities. There are many 
types of classical noise (thermal, ...) but, in addition, there exits intrinsic quantum noise. 
A wide variety of noise problems can be described in a universal manner by means of a 
generalized quantum Langevin equation (QLE) [lj]. This equation goes beyond the classi- 
cal Langevin equation |2| by incorporating a potential, both quantum and non-Markovian 
effects, the quantum fluctuation-dissipation theorem, |3| as well as the presence of a time- 
dependent external force (thus allowing a description of the evolution of an irreversible 
system). Our purpose here is to outline how one obtains the QLE by starting with a very 
general microscopic model. Next, we consider a variety of applications. Before proceed- 
ing, we will first give a brief history of different approaches to fluctuation and dissipative 
processes. 

As distinct from a dynamical system, a heat bath (reservoir) describes a subsystem with 
an infinite number of degrees of freedom which enable the system (i.e., the particle interacting 
with the heat bath) to relax in the course of time to a unique equilibrium state. A commonly 
used model of a heat bath is an infinite number of oscillators and the freedom of choice of 
the relevant masses and frequencies leads to a surprising diversity in the physical systems 
described, examples being blackbody radiation (BBR) and phonon heat baths and impurities 
in a metal (in the limit of the oscillator frequencies going to zero). 

The beginning of the subject is generally regarded as having started with the observations 
of Robert Brown, a Scottish botanist who observed the random motion of pollen grains 
immersed in a fluid |4]. No external forces were present and the temperature T was room 
temperature. Later work showed that the very irregular motion occurred in any suspension 
of tiny particles in a liquid medium. The explanation of these phenomena was given in a 
series ofpapers by Einstein [5j from 1905 to 1908. Coupled with the experimental work of 
Perrin £| this served to establish the atomic theory of matter since the irregular motion 
was clearly identified as being due to collisions with the molecules in the liquid. The term 
" Brownian motion" is now used in a generic sense to denote random motion and it covers a 
wide spectrum of phenomena from the motion of very fine particles suspended in gas to the 
motion of electrons in a BBR heat bath. Einstein's explanation of Brownian motion used 
a discrete time approach. In particular, his result that the diffusion constant D (which is 
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defined as one-half of the rate of change with time t of the mean-square displacement in the 
limit of large t is proportional to T/7, where T is the temperature and 7 _1 is the collision 
rate, is the first example of a fluctuation - dissipation (FD) theorem. 

Shortly after the work of Einstein, Langevin Q presented an entirely new approach to 
the subject which, in the words of Chandrasekhar pj, constitutes the "modern" approach 
to this and other such problems. The essence of Langevin's approach is a continuous time 
approach implemented by adoption of a stochastic differential equation, i.e., an equation 
for quantities which are random in nature. In other words, Langevin provided an elegant 
solution to the problem of generalizing a dynamical equation to a probabilistic equation. This 
was to be the start of a major new field of study with widespread implications in physics, 
chemistry, biology, and many other fields. The approach of Langevin was phenomenological 
in nature but its essential correctness has been verified by various microscopic studies. An 
essential feature of his approach was to separate the total force acting on a particle due to 
its environment into two parts: a frictional force and a fluctuation (random) force. These 
terms are very different in nature: The fluctuation term is basically microscopic in nature 
and has a time scale determined by the mean time between collisions whereas the time scale 
of the frictional force is proportional to the self-diffusion constant and is much larger. 

Another example of a fluctuation-dissipation theorem was provided by Nyquist [8] who 
showed that the random fluctuations in voltage and current across a resistor measured by 
Johnson 9] are determined by its impedance (the famous Johnson-Nyquist noise in electrical 
circuits). All of the aforementioned work was classical in nature, but in 1951, Callen and 

n 

Welton [3J presented their celebrated work on a quantum formulation of the FD theorem. 
Since a major shortcoming of the Langevin equation is its phenomenological nature, it is 
clearly desirable to obtain such an equation from microscopic considerations. Such a the- 
ory was given by Zwanzig [10] who considered the particle of interest to interact with an 
environment consisting of an infinite number of oscillators and obtained the usual classical 
high-temperature Brownian motion results. Turning to microscopic quantum dissipative 
phenomena, pioneering approaches to developing a quantum L ang evin equation for partic- 
ular problems appear in L article, of Se, ut : ky H a nd Lax jj. However, these a rtid e S 
were based on a Markovian approximation and, unfortunately, the work of Senitzky was 
marred by serious errors, as pointed out by Li et al.Q] Perhaps the most influential of 
the earlier articles is perhaps the oft-quoted article of Ford et al. ^ ] since it was the first 



3 



article in which the correct formulation of the quantum Langevin equation was indicated. 
Other important contributions include the work of Mori [15] , as well as Benguria and Kac 
jil and Ford and Kac leading up to the work of Ford, Lewis, and O'Connell (FLO) 
l| who, in an article entitled " Quantum Langevin Equation" , gave a detailed discussion of 
the problem. In particular, these authors presented the general form of the QLE consis- 
tent with fundamental physical requirements, in particular, causality and the second law of 
thermodynamics. Other approaches to dissipative problems include Kubo's linear response 
theory Q| , and the Feynman- Vernon use of path integrals which has been condider ably 
extended by Leggett [2^ and co-workers as well as many German investigators j^.j^ 
among others. However, the QLE approach is our method of choice since not only is it at 
least on a par technically with all the other approaches but we also consider it to be more 
physically appealing and simpler to execute. 

Sec. 2 is devoted to fundamentals. In particular, we show how the QLE can be obtained 
from a very general and basic Hamiltonian after which we proceed to show how this leads 
to results for important quantum commutators and correlations which form the basis for 
the derivation of observable quantities. The QLE derived in Ref. 1 pertains to a stationary 
process, in the sense that correlations, probability distributions, etc. for the dynamical 
variable x are invariant under time-translation (t — > t + to). It is this equation which is the 
basis of most of the applications which we discuss below. There is also a QLE for the initial 
value problem jl^] which has been used to obtain the most general master with explicit 
expression for the associated time-dependent coefficients, leading in turn to an explicit exact 
solution [2^. Sec. 3 discusses the Ohmic heat bath model and various applications and 
similarly for the radiation heat bath model in Sec. 4. Driven Systems are discussed in Sec. 5. 



II. GENERALIZED QUANTUM LANGEVIN EQUATION 

The QLE is a Heisenberg equation of motion for the coordinate operator x of a quantum 
particle of mass m moving in a one-dimensional potential V(x) and linearly coupled to a 
passive heat bath at temperature T. Whereas it has a very general form, as we pointed out in 
Ref. 1, it can be realized with a simple and convenient model, viz., the independent-oscillator 
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(10) model. The Hamiltonian of the 10 system is 

« = £ + + £ + im^to - - ./<«). (2.1) 

Here m is the mass of the quantum particle while rrij and Uj refer to the mass and frequency 
of heat-bath oscillator j. In addition, x and p are the coordinate and momentum operators 
for the quantum particle and qj and pj are the corresponding quantities for the heat-bath 
oscillators. Also f(t) is a c- number external force. The infinity of choices for the rrij and Uj 
give this model its great generality. In particular, it can describe nonrelativistic quantum 
electrodynamics (QED), the Schwabl-Thirring model, the Ford-Kac-Mazur model, and the 
Lamb model 

Use of the Heisenberg equations of motion leads to the QLE 

mx + f dt'n{t-t')x{t') + V\x) = F(t) + f(t), (2.2) 

J — oo 

where V'(x) = dV(x)/dx is the negative of the time-independent external force and fj,(t) is 
the so-called memory function. F(t) is the random (fluctuation or noise) operator force with 
mean (F(t)) = and f(t) is a c-number external force (due to an electric field, for instance). 
In addition, it should be strongly emphasized that "- the description is more general than 
the language — " [l( in that x(t) can be a generalized displacement operator (such as the 
phase difference of the superconducting wave function across a Josephson junction). 

Thus, the coupling with the heat bath is described by two terms: an operator-valued 
random force F(t) with mean zero, and a mean force characterized by a memory function 
fi(t). Explicitly, 

fx(t) = rrijU? cos(ujjt)9(t) , (2.3) 
j 

, with 6 {t) the Heaviside step function. Also 

F(t) = J2m j ^(t), (2.4) 

j 

is a fluctuating operator force with mean (F(t)) = 0, where q h (t) denotes the general solution 
of the homogeneous equation for the heat-bath oscillators (corresponding to no interaction). 
These results were used to obtain the results for the (symmetric) autocorrelation and com- 
mutator of F(t), viz., 

C FF {t-t>) = \(F{t)F(t') + F{t')F{t)) 
1 f°° 

= - / duRe[fl(uj + iO + )]hujcoth(huj/2kT)cos\uj(t-t')}, (2.5) 
n Jo ' 
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\F(t),F(t')} = — / duRe{jl(u + iO + )}u sinu(t - t'). (2.6) 
in Jo 

We note that both of these results are independent of V(x) and f(t). Here fl(z) is the 
Fourier transform of the memory function: 



Mz) = / dtfi(t)e izt . (2.7) 
Jo 

Kubo refers to (2.5) as the second fluctuation-dissipation theorem and we note 

that it can be written down explicitly once the QLE is obtained. Also, its evaluation requires 
only knowledge of Kefi(uj). On the other hand, the first fluctuation-dissipation theorem is an 
equation involving the autocorrelation oix(t) and its explicit evaluation requires a knowledge 
of the generalized susceptibility a(u) (to be defined below) which is equivalent to knowing 
the solution to the QLE and also requires knowledge of both Rejl(u) and lmp,(u). This 
solution is readily obtained when V(x) = 0, corresponding to the original Brownian motion 



problem [2j. As shown in Ref. 2^, a solution is also possible in the case of an oscillator. 
Taking V(x) = \Kx 2 = \miJ^x 2 these authors obtained the solution of the Langevin 
equation (2.2) in the form 

x(t)= f dt'G(t-t'){F(t') + f(t% (2.8) 

J — oo 

where G(t), the Green function, is given by 

1 f°° 

G(t) = — dtoa(uj + iO + )e~ luJt , (2.9) 
with a(z) the familiar response function 

^ ^ —mz 2 — iz[i{z) + K ^ 
Also, taking the Fourier transform of (J2.8)) . we obtain 

x(u)=a(u){F(u) + f(u)}, (2.11) 

where the superposed tilde is used to denote the Fourier transform. Thus, x(u) is the Fourier 
transform of the operator x(t): 

POO 

x(u) = / dtx(t)e lujt . (2.12) 
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It is also useful to note that the commutator, which is temperature independent, is given by 
the formula 2^ 

2ih r°° 

[x(t 1 ),x(t 1 + t)] = / dculm{a(uj + i0 + )} sin cut. (2.13) 

n Jo 

Whereas the autocorrelation and commutator of F(t) are independent of V(x) and f(t), 
this is not so in the case of x(t), as is obvious from (2.8) and (2.10). However, since most 
of our applications pertain to the case f(t) = 0, we will confine ourselves for the moment 
to this case and discuss later generalizations which are necessary when f(t) ^ 0. Thus, in 
particular, when f(t) = 

C(t - = X - (x(t)x(t') + x(t')x(t)) 

h r°° „ , _ , Uu 



- \ dujlm{a(u + iO + )}coth—— cos to(t-t f ). (2.14) 
7r jo 2kT 

lvof the 

mm 



We note that it is a function only of the time-difference t — t' . 

This is referred to by Kubo (l^Ql as the fluctuation-dissipation theorem of the first 



kind. Next, from (2.10), we see that 

Im a{uj) = uj\a{u))\ 2 Re jl{uj). (2-15) 

Thus, (2.14) may be expressed in the form 

h r°° 



C{t) 
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f°° ( Tiuj \ 

/ doj u;|a;(c<;)| 2 Re coth cos tut. (2.16) 

Jo \2kT J 



The mean square displacement of the quantum particle in a dissipative environment, s(t) 
say, plays a key role in many of our subsequent discussions. In particular, it determines the 
diffusion time through a system of interest. Thus using (|2.14|) . we obtain 



s(t) = ([x s (t)-a: s (0)] 2 ) 



= 2{C(0)-C(t)} 

2h r°° htu 
= — dulm{a(uj + i0 + )}coth— — (1 -cosut). (2.17) 
7r jo 2kT 

For specific applications, the key ingredient to be identified is the nature of the heat bath 
which is simply characterized by Re jl(uj + i0 + ), the spectral distribution of the memory 
fuction, which is a positive real function and given explicitly by the relation 

Re \fi(u + i0 + )] = ^Y!,mj(J$[5(u - u,) + 5(to + Uj)). (2.18) 

j 



Next, we summarize how the above tools have been used to obtain thermodynamic quanti- 
ties, especially results which have been obtained only by the QLE method. The key quantity 



calculated 



24j is the free energy F, which is a thermodynamic potential from which other 



thermodynamic functions can be obtained by differentiation. 

The system of an oscillator coupled to a heat bath in thermal equilibrium at temperature 
T has a well-defined free energy. The free energy ascribed to the oscillator, F(T), is given 
by the free energy of the system minus the free energy of the heat bath in the absence of 



the oscillator. This calculation was carried out by two different methods 2^ |. 
the " remarkable formula" : 



1 f°° 

F(T) = - / duf(u,T) x Im 

7T JO 



d\oga(uj + i0 + ) 



leading to 



(2.19) 



where f(u,T) is the free energy of a single oscillator of frequency u, given by 

I —huj 



f(u,T) = KT log 



exp 



kT 



(2.20) 



Here the zero-point contribution (4^) has been omitted. In particular, this result was used 
to obtain the energy and the entrophy ascribed to the osillator. 

Finally, we note that magnetic field effects on all of the above quantities have been 
calculated j^,!^. In particular, it was shown that the effect of an arbitrary heat bath on 
Landau diamagnetism is always such as to reduce the magnitude of the magnetic moment 
without changing its diamagnetic character 

Thus, we now have all the tools necessary for applications. For each particular application, 
one must first decide on the correct specification of the relevant heat bath. The case of 
constant friction, the so-called Ohmic model, is of special interest since it is the simplest 
model but yet gives a good description of many physical systems. A generalization of this 



model, involving an additional parameter, is the single relaxation time model 



Also of 



interest is the Debye mode] 29J. The blackbody radiation (BBR) heat bath model 
is of fundamental interest since the relevant Hamiltonian H in this case is the universally 
accepted H of quantum electrodynamics. Thus, we refer to the BBR model as "the rosetta- 
stone of heat bath models" since it leads to readily-checked predictions and it is the basis of 
many of the applications which we consider below. Hence, in Sec. 3, we consider the Ohmic 
model and various applications and we do likewise for the BBR model in Sec. 4. 
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III. OHMIC HEAT BATH MODEL AND APPLICATIONS 



The Ohmic model is defined by the choice 

fi(t)=(5(t), t>0, (3.1) 

so that 

//(a?) = C = m l — constant. (3.2) 

Also, it follows that 

m{x + ^x + ul) = F{t). (3.3) 

This corresponds to the original form of the Langevin equation except that now the classical 
quantities x(t) and F(t) are operators and we have also included an oscillator potential. 
Since the past motion does not appear, one says there is no memeory. On the other hand, 
the quantum- mechanical process is not Markovian since Cppif — £Q is n °t proportional to 
Cpp(t — t') a 5(t — t') except in the classical limit where h — > [lj. Using (3.3) in (2.10) 
leads to the result 

a(u) = {m[-uj 2 + uj 2 -iuj]}' 1 , (3.4) 
from which it follows, using (2.9), that the green fucition is given by 

G{t) = e-& ^ (3.5) 



where 



Also, from (2.16), we obtain 



^i = H 2 -(|) 2 }^- (3-6) 



tvy_ r cucoth(^) 

271W 

We now consider various applications 



c (t) = i 2 "^;;;t 2 2 - ( 3 - 7 ) 



A. Brownian Motion 

This is the original problem of diffusion of a free particle through a medium characterized 
by a dissipative decay rate 7. Since here u = 0, the results given above reduce to 
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a(u) = {m[-uo 2 -iuj]}-\ (3.8) 

<m = (3.9) 

£7( t ) = «3L f rot ^> e"«. (3.10) 



27rm j-oo cj(u; 2 + 7 2 ) 
Thus, (2.17) gives 

h'j f°° , (1 - cos cut) / hu \ , ^ 

s(t =— / ^ 1 7 2 ^ 2 ( ooth — . 3.11 

In the classical high temperature limit (corresponding to the Brown experiment), 
coth (#)— (ig^so that 

s(t) = r |e" 7i - 1 + 7 t| , kT»h~/. (3.12) 

For long times 

2A;T 

s(t) -> t, 7* » 1, (3.13) 

m7 

which is the familiar Einstein relation from Brownian motion theory. However, for short 
times (which are more characteristic of decoherence decay times) 

kT 

s(t) -> — t 2 , 7t « 1, (3.14) 
m 

independent of 7. Thus, in this simple case, the familiar Einstein diffusion coefficient 

1 kT 

D = -m = —, (3-15) 

2 7777 

is obtained. 

In particular, the situation is very different for low temperatures (kT << Hj), in which 
case the main contribution is from the zero-point (T = 0) oscillations of the electromagnetic 
field. This calculation has recently been carried out |30( leading to the result 



s(t) = -^- 2 t 2 \log C - + lE - 3 -}, (3.16) 
where 7^ = 0.577215665 is Euler's constant. 
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B. Noise in gravitational wave detector suspension systems 



Crucial to the detection of gravitational waves is the mechanical system which is used to 
measure the relative displacements of suspended mirrors [interferometric experiments, such 
as the Laser Interferometric Gravitational Wave Observatory (LIGO) or the displacement 
of a resonant mass antenna (bar experiments)]. Here, we present a general framework for 
describing noise effects in detectors systems. Our results apply to a very general dissipative 
environment but for defmiteness we will concentrate mostly on the LIGO-like test mass 
suspensions and on bar detectors. 

The crucial quantity to be calculated is the ensemble average of the square of the dis- 
placement due to noise, (x 2 (t)), which is simply equal to C(0), as is clear from (2.14). Hence, 
from (2.16), 

h r°° ( hu \ c°° 

l x 2 (t)) = - / duu | a(u) | 2 Rejitu) coth — — = / P{u)du, (3.17) 
7r Jo \2kT J Jo 

where 

P{u) = 5« | a(u) | 2 Refi{u>) coth , (3.18) 

is the power spectrum of the coordinate fluctuations. 

In the case of resonant bar detectors, Ji{uj) is taken to be rwy, where 7 is a constant. 
Thus, in particular, if we take the weak-coupling limit ii) ■€ 7 and substitute (3.4) into 
(3.17), we immediately find the well-known high-temperature result (x 2 (t)) = and the 
zero-temperature result (x 2 (t)) = ( ) [3JJ. However, in the case of nonresonant LIGO 
detectors, which are responsive to a range of frequencies, the frequency dependence of Jl{ui) 
is essential. In practice, because of the complexity of the detector system, the practical 
procedure will be to fit the experimental results by some Jx{u). This has the advantage 



reat bath. 



that we know a lot about the properties of ju(u;), regardless of the nature of the 
In particular, it is independent of the external potential and the temperature [l|. Further 
details may be found in fl. 

C. Resistively shunted Josephson junctions 

Our purpose here is to show how some of the principal well-known theoretical results 
follow simply from the quantum Langevin approach and also how they may be generalized. 
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This will also provide an example of our statement that "- the operator x is the quantum 
Langevin equation - can be a generalized displacement operator -". For an ideal junction 



the current is given by the Josephson equation [321 ]: 



I = losing, (3.19) 

where <fi is the phase difference of the superconducting wave function across the junction 
and Ic is the critical current. The voltage across the junction is 

V = ^J>, (3.20) 

where ^ is the quantum of flux. A real junction can be viewed as a capacitance C and a 
shunt resistance R in parallel with an ideal junction. The current is then the sum of the 
ideal junction current, given by (3.19), the current through the capacitor, Q = CV, and the 
current through the resistor, / = ~. The junction voltage is still given by (3.20). The basic 
equation of motion of the junction can therefore be written 

This is of the form of a quantum Langevin equation with mass and friction constant 



and with potential 



U{<f>) = ~{I<f> + I c co&4>). (3.23) 
2e 



Thus, analogous to (3.17), the ensemble average of the square of the phase due to noise is 

<r>= ^r^ coth (^) % _ (3.24) 
7rChJo \2kT ) (Uq — u 2 ) z + u z Y 

where 

1=W- ^W^o-P)*. (3.25) 
In the limit of large shunt resistance (weak coupling limit), 7 <C ujq, this becomes 

This weak coupling limit corresponds to the expression for the phase fluctuations obtained 
by Josephson |33j • The power spectrum of the voltage fluctuations is readily obtained from 
(3.24) using (3.20). Further details may be found in (34). 
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D. Environmental Effects on Nanosy stems 



The study of nanosystems is now burgeoning and, because of the small dimensions in- 
volved, environmental effects play an important role. To illustrate, we will focus our discus- 
sion on the analysis of charge fluctuations on small-capacitance tunnel junctions. Follow- 
ing Ingold and Nazarov 3^|, we consider a junction of capacitance C carrying the charge 
Q = CV where V is the voltage across the junction. The external circuit (the environment) 
is described by its impedance Z(u>) = 7M, where is the current. The phase differ- 
ence across the junction is obtained from (3.20) [except the 2e — > e since here we are not 
considering a superconducting tunnel junction] to get 

0(t) = £ f dt'V{t r ). (3.27) 

Turning to a quantum picture, it follows that (- )0 and Q are the electrical equivalent of the 
mechanical quantities x and p (see table I of Ref. 35 for more details). It readily follows 
that there is a fluctuation-dissipation theorem describing the relation between the charge 
fluctuation q(t) and the dissipation described by Z{uj). Explicitly, we have 

TiujC 2 , hu 



< q (t) >= / du coth' 

J 7T 



-Re 



iujC + Z~ 1 (uo) 



(3.28) 



2k B T 

This result is applicable to any kind of external circuit attached to a small junction. In 
particular, this result was used to calculate the effect of quantum smearing of Coulomb 

was shown that, in the weak coupling limit, 
the QLE theory is related to the well-known generalized Landauer formula j38|. Finally, 
we note that big efforts are currently being made toward developing single charge transfer 
devices and, in that context, we have used (3.27) to calculate environmental efforts on a 
single electron box [39]. 



E. Quantum Transport for a many-body System 

The problem of electrical conductivity in solids has been studied by a large number of 
theoretical methods. For the most part, they are either based on the Boltzmann transport 
equation for the distribution of electrons, or they start with the Kubo formula in linear- 
response theory. In spite of their successes, those two methods are known to be impractical 
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tools for calculating the electrical transport properties of systems with high-order impurity 
scattering, among other disadvantages. Thus, we were motivated to treat the quantum 
transport of an interacting system of electrons, impurities, and phonons, in a time- dependent 
electric field, by using the QLE in which the system is shown to be equivalent to a quantum 
particle in a heat bath [40]. The center of mass of the electrons acts like a quantum particle, 
while the relative electrons and phonons play the role of a heat bath. They are coupled 
through the electron-impurity and electron-phonon interactions. After eliminating the heat- 
bath variables, the equation of motion for the quantum particle is written in a form of 
a QLE, with a memory term which reflects the retarded effects of the heat bath on the 
quantum particle. The evaluation of the memory term immediately leads to a result for the 
susceptibility from which we can calculate the conductivity directly, in contrast to Kubo- 
type calculations which require the evaluation of correlation functions as an intermediate 
step. Our results were then applied successfully to transport in 3, 2, 1, and dimensional 
systems and we refer to our review j^lj for details. As emphasized in j4l|, a big advantage 
of the QLE approach is that it gives a natural separation between the conductivity and the 
noise. 

F. Unruh radiation 

It is generally accepted that a system that undergoes uniform acceleration with respect 
to the vacuum of flat-space-time will thermalize at a temperature (the so-called Unruh 



temperature) that is proportional to the acceleration 42] . However, the question of whether 
or not the system actually radiates is highly controversial. Thus, we were motivated to 
carry out an exact analysis of the problem using a generalized quantum Langevin equation to 
describe the motion of an oscillator (the detector) moving under a constant force and coupled 
to a one- dimensional scalar field (scalar electrodynamics). We conclude that the system does 
not radiate despite the fact that it does in fact thermalize at the Unruh temperature Q|. 
The essence of the calculation was to show that the equation of motion of the oscillator 
has the same form as that given in (3.3) except that now the autocorrelation given in (2.5) 
not only corresponds to an Ohmic heat bath but the relevant temperature is the Unruh 
temperature which, of course, is proportional to the acceleration. However, due to the fact 
that there is no term analogous to the external force f(t) appearing (2.2), it immediately 
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follows that there is no radiation, as confirmed by a detailed calculation. What we found is 
that the uniform acceleration does not give rise to a constant force in the equation of motion 
but, instead, it only appears explicitly in the autocorrelation of the fluctuation force. In 
other words, in this case, the constant external acceleration does not cause radiation but 
simply accentuates the noise. 



IV. BLACKBODY RADIATION HEAT BATH MODEL AND APPLICATIONS 

The relevant Hamilton H is the universally accepted H of quantum electrodynamics 
(QED) from which it follows Q,Q that 

Re[Ku + i + )} = ^fl (4.1) 

where the quantity fk is the electron form factor (Fourier transform of the electron charge 
distribution). In other words, we have allowing the electron to have structure. 

The physically significant results for this model should not depend upon details of the 
electron form factor, subject, of course, to the condition that is be unity up to some large 
frequency Q and falls to zero thereafter. A convenient form which satisfies this condition is 

fi = OT' (4 ' 2) 

Using this in (4.1), the Stieltjes inversion formula gives 

We see here a manifestation of the general feature that the memory function is independent 
of the external potential and the particle mass. 

As essential aspect of QED theory is the necessity for mass renormalization. Thus the m 
occurring in the QLE is actually the bare mass and the renormalized (observed) mass M is 
given in terms of the bare mass m by the relation 

M= m + 2 f Q =m + T e QM, (4.4) 

where 

2e 2 

6xl0~ 24 sec. (4.5) 



3Mc 3 
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Next, taking the inverse Fourier transform of (4.3) leads to the result 

fi(t) = Mn 2 r e [26(t) - fiexp(-fit)], (4.6) 

which makes manifest the non-Markovian nature of the motion. We now consider various 
applications. 

A. Elimination of Runaway Solutions in Electron Radiation Theory 

The fact that ji(uj) is a positive function (which follows from the second law of thermody- 
namics) is equivalent to the demand that all the poles of a(u>), the generalized susceptibility, 
must lie in the lower half of the complex plane (which also follows from the principle of 
causality). This implies that m > and that 

Q, < r- 1 = 1.6 x 10 23 s -1 . (4.7) 

If one selects for Q its largest permissible value: namely r^ 1 ( corresponding to letting m — > 0, 
which is also equivalent to choosing the closest approach to a point electron consistent with 
causality), then one obtains the classical equation of motion 

Mx = f{t)+r e f{t). (4.8) 

This result has a variety of desirable properties not exhibited by the Abraham-Lorentz 
equation. Generalizations which include quantum effects and the presence of a p otential 
may be found in [^jQ]- The Larmor formula also requires generalization |46| . . A 
relativistic generalization of (4.8) is given in [48] along with a proof that a constant force 
does not cause radiation. 



B. Rydberg Atom Level Shifts due to Blackbody Radiation 



Atoms coupled to a bath in thermal equilibrium at temperature T are best described in 
terms of the free energy 1^1 and we have argued that this is the quantity actually 

measured by experiment [51|. Thus, using (2.19), we found that, in the high temperature 
limit, where kT is large compared with the level spacing, the shift in free energy is given by 

«w (49) 

in agreement with the experimental results jsil ]. 
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V. DECOHERENCE 



The superposition principle (and related work on Schrodinger cats), entanglement, and 
the quantum-classical interface are at the cutting-edge of topical research [52j. Since su- 
perposition states are very sensitive to decoherence, reservoir theory has attracted much 
recent interest. Recognizing that conventional master equation approaches are often not ad- 
equate, some investigators have used path integral methods. However, we have found that 
the simplest and most physically appealing approach to the problem is via use of the QLE, 
supplemented by use of the Wigner distribution function (WDF) for the study of a system 
with an infinite spectrum of states and by use of the spin polarization vector in the case 
of spin systems |54j - For the investigation of decoherence phenomena, we have found that 
it is important to distinguish between two different physical scenarios: namely (a) complete 
entanglement between the quantum particle and the heat bath at all times 55j and (b) the 
system in a state in which the oscillator is not coupled to the bath at, say, t — and such 
that the bath is in equilibrium at temperature T [5j|. Thus, it takes a characteristic time 
of the order of 7" 1 (where 7 is a typical dissipative decay time) for the complete coupling 
to occur and for the whole system to come into thermal equilibrium. It is notewor thy that 
"decoherence without dissipation "can occur [55-57] but not for low temperatures 58]. 

Scenario (a), the "entanglement at all times" calculation utilized quantum probability 
distributions (which are related to Wigner distributions) in conjunction with results obtained 
by use of the stationery solution to the QLE, given in (2.2). On the other hand, for the 
scenario (b) case, it is necessary to take into account initial conditions. The Langevin 
equation for the oscillator with given initial values is given by 

mot + /* dt'fj,(t - t')x(t') + Kx = -u(t)x(0) + F(t), (5.1) 
Jo 

and the general solution is given by 

x(t) = mG{t)x{0) + mG(t)x{0) + X(t), (5.2) 

where we have introduced the fluctuating position operator, 

X(t) = f dt'G(t-t')F(t'). (5.3) 
Jo 

If we assume that at t = the system is in a state in which the oscillator is not coupled to 
the bath and that the bath is in equilibrium at temperature T, we find that the correlation 
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and commutator are the same as those for the stationary equation. It was then possible to 
show 5(| that one could write the Langevin equation (|5.1|) in the form of an equation that 
is local in time with time-dependent coefficients: 

x + 2T(t)x + n 2 (t)x = — F(t), (5.4) 

m 

where explicit expressions for T(t) and Q(t) were obtained in terms of G(t). Furthermore, 
it was shown that these results constitute in essence a derivation of the HPZ exact master 
equation [53 with explicit expressions for the time-dependent coefficients. It is also notable 
that in this scenario one needs to assume that the initial temperature of the particle is the 
same as that of the heat bath in order to obtain "decoherence without dissipation". 

VI. DRIVEN SYSTEMS 

As remarked in Sec. II, the position autocorrelation, as distinct from the fluctuation force 
autocorrelation, is modified by a presence of an external force f{t), as discussed in detail in 
It is convenient to now write 

x(t) = x s (t) + x d , (6.1) 

and 

C(t) = C {t)+C d (t), (6.2) 

where x d is the "driven" contribution due to the external force f(t) and x s is the contribution 
due to the fluctuation force F(t). Here we have introduced a subscript s to emphasize that 
x s (t) is a stationary operator-process, in the sense that correlations, probability distributions, 
etc. for this dynamical variable are invariant under time-translation (t — > t + to). In 
particular, the correlation Co(t) is given bv the right-side of (2.14) and is a function only of 
the time-difference t — t'. Furthermore, 6J| 

C d (t,t') = {x(t)){x(t')), (6.3) 

where (x(t)) is the steady mean of the driven motion. Based on these results, we extended 
the calculation of the well-known Burshtein-Mollow spectrum of resonance fluorescence to 
the case of non-zero temperature. In the high temperature limit, nT » hu , where Hlvq 
is the resonant energy for the two- level (atom), we found that the decay rate is increased 
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by a factor Our calculation was based on use of the Lax formula for two-time 

correlations but, as we also pointed out, this formula is applicable only for weak coupling 



and for frequencies near a resonance frequency since, more genera 
regression theorem cannot be generalized to the quantum domain 



ly, t he Onsager classical 
6z 



A further application of these results was the consideration of decoherence phenomena in 
the presence of an external field jf^J, 6^|. This work was motivated by the recent interest 



in engineered reservoirs |65l |. In particular, Myatt et al. [66] used a linear Paul trap to 
confine single Be ions in a harmonic potential and then prepared various superposition 
states. Next, they induced decoherence by coupling the single ion to a reservoir which 
they controlled in various ways. Such a reservoir gives rise to an external force f(t) in the 
equation of motion of the system. The calculations in j^J, (m| led to the conclusion that a 
non-random external field does not give rise to decoherence whereas, by contrast, a random 



field does. In particular, the experiments of Myatt et al. 66] used a 5-correlated force. 
The existing experiments verified the familiar result that the decoherence decay time is 
inversely proportional to the square of the separation d of the superposition components. 
However, this is a familiar result predicted by the plethora of papers dealing with the f(t) = 
situation but it does not give information on the dependence of on the parameters of the 
externally-superimposed reservoir. Thus, we will have to await come experimental data in 
order to compare with existing theory. 
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